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Abstract. We present a detailed spatial and dynamical analysis of the central ~ 2.2 Mpc region of the 

galaxy cluster Abell 521 (z=0.247), based on 238 spectra (of which 191 new measurements) obtained at the 3.6 m 
Telescope of the European Southern Observatory and at the Canada-France-Hawaii Telescope. From the analysis 
of the 125 galaxies that are confirmed members of the cluster, we derive a location (“mean” velocity) of Cm = 
74019j:“^ km/s and detect a complex velocity distribution with high velocity scale (“dispersion”, Sbi = 13251jqq 
km/s), but clear departure from a single Gaussian component. When excluding a possible background group of 
four galaxies, the velocity dispersion remains still large (~ 1200 km/s). The general structure of the cluster follows 
a North-West/South-East direction, crossed by a perpendicular high density “ridge” of galaxies in the core region. 
The Northern region of the cluster is characterized by a lower velocity dispersion as compared to the whole cluster 
value; it hosts the BCG and a dynamically bound complex of galaxies, and it is associated with a group detected 
in X-ray (Arnaud et al. 2000). This region could be in a stage of pre-merger onto the main cluster. The small 
offset (~ -I- 250 km/s) in the mean velocity of the northern region as compared to the whole cluster suggests that 
the merging occurs partly in the plane of the sky. These results, taken together with the fact that most of the 
clumps detected on the isodensity maps, as well as the early-type galaxies and the brightest ones (L> L* *) are 
aligned, suggest that this North-West/South-East direction is the preferred one for the formation of this cluster. 
The central high dense region (“ridge”) shows a lower velocity location (Cbi = 73625 l/ 3 gQ km/s) and significantly 
higher scale (1780jlj42 km/s) as compared to the whole cluster values. This is due to the presence of a low-velocity 
group of galaxies with a high fraction of emission line objects. This can be explained in a scenario in which a 
merging of subclusters has recently occurred along the direction of the “ridge” with a significant component along 
the line of sight. The low-velocity group would then be a high-speed remnant of the collision which would have 
also triggered an episode of intense star formation responsible for the large fraction of late-type objects in this 
region. 
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1. Introduction 

In the hierarchical model of structure formation, galaxy 
clusters are supposed to form by merging of units of 
smaller mass. Analysis of statistical samples of galaxy 
clusters have shown that a high percentage of clusters 
with substructures is detected even at low redshift, im¬ 
plying that clusters are still today undergoing the process 
of formation (Geller & Beers 1982, Dressier & Shectman 
1988, Jones & Forman 1992). Moreover, quantifying pre¬ 
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cisely the amount of morphologically complex clusters al¬ 
lows one in principle to constrain directly the cosmologi¬ 
cal model through the density parameter 0^ (Richstone 
et al. 1992, Mohr et al. 1995). This analysis is however 
hampered by the existing uncertainty in the rate at which 
substructure is erased (Kauffmann & White 1993, Lacey 
& Cole 1993). The second difficulty is that sub-clustering 
can affect the various quantities observable, such as the 
projected distribution of the galaxies and of the gas, the 
velocity distribution of the galaxies and the temperature 
structure of the gas, not necessarily at the same level, 
leading to sometimes different conclusions. 

Detailed studies of individual complex galaxy clus¬ 
ters at different wavelengths is a complementary analysis 
which allows one to obtain details of the scenario of forma¬ 
tion of these objects, and the physical processes necessary 
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to explain the observed distribution (Flores et al. 2000, 
Donnelly et al. 2001, Mohr et al. 1996, Bardelli et al. 1998, 
Rose et al. 2002, Berrington et al. 2002, Czoske et al. 2002, 
Valtchanov et al. 2002). In this paper, we will concentrate 
on the dynamical analysis of the merging cluster Abell 521 
which has been targeted for its outstanding properties. 

Abell 521 is a rich (R=l) Abell cluster, first detected 
in X-ray with HEAOl (Johnson et al. 1983, Kowalski et 
al. 1984). It was also suspected (Ulmer et al. 1985) to 
form a binary cluster together with its nearest neighbor 
on the sky, A518, but no clear evidence for gas interaction 
was found between the two clusters. Radio observations 
in the region of this cluster (Hanisch et al. 1985) have 
also shown a high fraction of radio sources with projected 
distances to the center compatible with these objects be¬ 
ing cluster members. More recent data, both in X-ray and 
optical (Arnaud et al. 2000, Maurogordato et al. 2000) 
provided a more detailed analysis of the properties of the 
galaxies and gas distributions in this cluster. Imaging in 
X-ray (ROSAT/HRI) has shown a gas morphology with 
two peaks which can respectively be associated with a 
diffuse main cluster, and a compact less massive group 
in the northern region, suspected to be in pre-merger 
stage with the main cluster. The projected galaxy den¬ 
sity distribution in the central 2.2 hya”^ Mpc has a very 
anisotropic morphology, as it exhibits two high density 
filaments crossing in an X-shape structure at the barycen- 
tre of the cluster. A severe gas/galaxy segregation stands 
out. The brightest cluster galaxy is offset of the cluster 
barycentre, and lies in the region of the X-ray northern 
group. Multi-object spectroscopy at ESO/EFOSC2 and 
CFHT/MOS led to the determination of the mean red- 
shift of the cluster, z=0.247, and of its velocity dispersion 
Sbi = 1386 km/s as measured from 41 members. However, 
this very high value of the velocity dispersion could be af¬ 
fected by the presence of substructures. Its value is also 
high compared to the temperature of the X-ray gas mea¬ 
sured with ASCA (T=6.3 KeV, Arnaud et al. 2000). These 
results imply that this cluster is undergoing strong dynam¬ 
ical evolution. This motivated new observations in order 
to better characterize the merging scenario, in particular 
through additional multi-object spectroscopy. A project 
of wide-held multicolor imaging in hve bands is under 
progress, and will address the large scale environment of 
the cluster using photometric redshifts (Ferrari et al. in 
prep.). In this paper, we analyze in detail the velocity dis¬ 
tribution within the central 2.2 hya”^ Mpc of the cluster 
with new data obtained at the ESO 3.6 m telescope (191 
new redshifts measured). Section 2 briehy describes the 
observations and the data processing techniques, as well 
as the level of completeness achieved. In Sect. 3, we per¬ 
form a general analysis of the ID velocity distribution, 
test for departures from an unimodal gaussian, and ht a 
partition in three velocity groups. In Sect. 4, correlations 
between the structures identihed in velocity space and in 
projected coordinates have been looked for. Variations of 
dynamical properties with absolute luminosity, color, and 
spectral type have been addressed in Sect. 5. In Sect. 


6, our new data are used to elaborate the more plausi¬ 
ble scenario of occurrence of the various merging events 
within this particularly complex cluster. All numbers are 
expressed as a function of h 75 , the Hubble constant in 
units of 75 km/s/Mpc. We have used the ACDM model 
with = 0.3 and Da = 0.7, then 1 arcmin corresponds 
to ~0.217 h 75 “^ Mpc in the following. 

2. The data 

2.1. Observations and data reduction 

New data have been obtained through a campaign of 
multi-object spectroscopy at the ESO 3.6 m telescope 
(3 nights in October 1999, and 2.5 nights in December 
2000). We used the ESO Faint Object Spectrograph 
and Camera (EFOSC2) with grism#04, whose grating of 
360 line mm“^ leads to a dispersion of 1.68 A/pixel, and a 
wavelength coverage ranging from 4085 to 7520 A. The de¬ 
tector used was the EFOSC2 CCD Loral/Lesser^40, with 
an image size of 2048 x 2048 (we made a 2 x 2 binning, 
in order to improve the signal to noise ratio), and a pixel 
size of 15 X 15 /rm. 

During the run of October 1999 we achieved a spectral 
resolution of FWHM^18.5 A, while in the second run, 
as a smaller punching head was available (1.35" instead 
of 1.8"), the resolution was improved to FWHM'^12.5 A. 
The total integrated exposure time was 9000 s for each 
frame, split in at least two exposures to eliminate cosmic 
rays. After each science exposure, a Helium-Argon lamp 
exposure was systematically taken for wavelength calibra¬ 
tion. 

Data have been reduced with IRAfQ, using our au¬ 
tomated package for multi-object spectroscopy based on 
the task “apall”. Radial velocities were determined us¬ 
ing the cross-correlation technique (Tonry & Davis, 1981) 
implemented in the RVSAO package (developed at the 
Smithsonian Astrophysical Observatory) with radial ve¬ 
locities standards obtained from the observations of late- 
type stars. We have obtained 191 new spectra. Among 
these, 29 are stars, while 109 are identified as galaxy spec¬ 
tra with a signal to noise ratio sufficient to obtain ra¬ 
dial velocity measurement with a parameter R of Tonry 
& Davis greater than 3. The remaining 53 objects have a 
poor velocity determination. 

We list in Table |I^ our new velocity measurements. 
The columns read as follows: Col. 1: identification number 
of each target galaxy; Col. 2: run of observations (October 
1999=ES01, December 2000=ES02); Cols. 3 and 4: right 
ascension and declination (J2000.0) of the target galaxy; 
Cols. 5 and 6: best estimate of the radial velocity and as¬ 
sociated error from the cross-correlation technique (those 
values have been set to “-2” if the object is a star and to 
“-1” if we have no redshift information); Col. 7: a quality 

^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper¬ 
ative agreement with the National Science Foundation 
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Fig. 1. Central field (20’xl0’) of the cluster. In black, galaxy iso-density contours for B<27 and I< 20; the lowest 
iso-density contour corresponds to Icr level above the mean density in the field, the contours are spaced by 0.5cr. 
Black circles correspond to the 125 galaxies of the cluster with spectral quality flag=l. The subgroups identified in 
the projected density map and the two main directions of the cluster (SI and S2) have been shown. 


flag for the redshift determination: l=good determination 
(R > 3), 2=uncertain determination, 3=very poor deter¬ 
mination, 4=failed spectra, Col. 8: a listing of detected 
emission lines. 

2.2. The spectroscopic sample 

In the following analysis, the new set of spectroscopic data 
presented above has been combined to our original sample 
(Maurogordato et al. 2000, 47 objects added to our new 
catalogue) resulting in a sample of 209 galaxies and 29 
stars. 

B- and I-band imaging (CFH12k) taken as part of our 
multicolor imaging program were used to build a color 
catalogue in the central region of the cluster surveyed 
by spectroscopy. This allowed to associate B- and I-band 
magnitudes for most galaxies of the spectroscopic sample 
(187), except for 22 objects located in the gaps between 
the chips of the camera or strongly blended. Fig. ^ dis¬ 
plays the central field of the cluster, where galaxies with 
secure redshift determination have been circled. The iso¬ 
density map of the projected distribution of galaxies with 
I-band magnitude I<20 is also displayed (derived using 
the Dressier algorithm; Dressier 1980). The preferential 
directions SI and S2 observed in Arnaud et al. 2000 are 
indicated. The density structures detected at more than 


5cr level are also indicated: six groups called gl to g6, a 
group around the BCG, and the so-called “ridge” struc¬ 
ture corresponding to SI. 

In Fig. H we have plotted the I-band magnitude dis¬ 
tributions of: the galaxies of our magnitude/velocity sam¬ 
ple (187), those with a very good redshift determination 
(141), and, among them, those belonging to the clus¬ 
ter (113). In Fig. ^ we show the ratio of the number 
of objects with measured velocities to the total number 
of galaxies detected within the central 10x10 arcmin^ 
of the field as a function of the I-band magnitude. We 
reach a general level of completeness for spectroscopy of 
50% at Iab = 19.5, which drops at 30% at Iab = 20.5. 
However, these values are strongly affected by several in¬ 
complete fields at the periphery of the cluster, as the cen¬ 
tral dense regions of the cluster have been much better 
sampled. In fact, we have divided the spectroscopic field 
in 2.5x2.5 arcmin^ cells (Fig. |^), and measured the de¬ 
gree of completeness in each cell for three different cuts 
in I magnitude (Iab = 19,20,21). The obtained values 
are shown in each cell of the corresponding isodensity 
maps. At Iab = 19, we have an excellent velocity sampling 
(^75-^80% completeness) in the North-West/South-East 
main structure of the cluster, which drops at ~60% at 
Iab = 20, and at ~40% at Iab =21. 
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Fig. 2. I-band magnitude distributions of the galax¬ 
ies of our spectroscopic sample (187 objects - solid 
line), of all the galaxies with good velocity determi¬ 
nation (141 - dot), and, among them, of those be¬ 
longing to A521 (113 - dash). 
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Fig. 3. Velocity completeness for different cuts in I- 
band magnitude in the central 10x10 arcmin^ field 
covered by spectroscopy. 


3. The velocity distribution 

3.1. Global analysis of the velocity distribution 

We have analyzed the general behavior of the velocity dis¬ 
tribution with the ROSTAT package (Beers et al. 1990). 
For this purpose, and in all the following analysis, we have 
used only the 125 objects with quality flag=l (secure red- 


1 1 1 1 1 1 1 1 1 

20.0% 

1 1 1 1 1 1 1 1 1 

55.6% 

1 1 1 1 1 1 1 1 1 

33.3% 

1 1 1 1 1 1 1 1 1 

28.6% - 




■ - 

66.7% 

77 8% 

86fc 

50.0% - 

: F 

yi- 






- 

-n 

: It' 

73.5% 

T 


■ 

0.0% - 

1 1 1 1 1 1 1 1 1 

(ji 

V 

30.0% 

0.0% 

0.0% - 



1 1 1 1 1 1 1 |i 1 

8 6% 

11 1 1 1 1 1 1 1 1 

18.2^ 

1 1 1 1 1 1 1 1 1 

24 0% 

1 1 1 1 1 1 1 1 1 

1 1 2% 




: 1 

P 

- 



! 


- 

*=■ 

i 

■ 

- 

. 

1 39.7%^ 

1 m^2% 

15.2% - 

35.4% 


20.6% 

0.0% - 


' ■■ 


- 

_ LB 1 1 





! 28.0% 

5.4% 

0.0% - 

. 

. 

. 

. 


Fig. 4. Completeness factor is superimposed to the 
isodensity maps in the field selected for spectroscopy 
(lO’xlO’). Each cell of the grid covers 2.5’x2.5’. 
Different magnitude cuts have been considered; from 
top to bottom: I<19, I<20, I<21. 
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Fig. 5. Velocity histogram within Abell 521 obtained 
from the 125 Q.F.=1 members of the cluster, with a 
binning of 500 km/s. A gaussian function with the 
velocity distribution derived with ROSTAT is super¬ 
imposed. 


shift) in our dataset; their velocity histogram is shown in 

Fig. |. 

In the case of large number of redshifts (100<n<200) 
as in our situation, the best choice to estimate location 
(“mean” velocity) and scale (velocity “dispersion”) is the 
biweight estimator (Beers et al. 1990), as it provides the 
best combination of resistance and efhciency across the 
possible contaminations of a simple gaussian distribution. 
We find a location Cbi = 740191^25 km/s and a scale 
Sbi = 132511^5 km/s. These results are in good agreement 
with those obtained from the analysis of our former sample 
of 41 galaxies (Maurogordato et al. 2000), and confirm the 
apparently high value of the velocity dispersion in this 
cluster. 

We have also analyzed the higher moments of the dis¬ 
tribution, in order to look for possible deviations from 
gaussianity that could provide important signature of dy¬ 
namical processes. For all the following tests the null 
hypothesis is that the velocity distribution is a single 
Gaussian. The traditionally used shape estimators are kur- 
tosis and skewness; in addition, we have computed the 
asymmetry and tail indices (Al and TI), which also mea¬ 
sure the shape of a distribution, but are based on the order 
statistics of the dataset instead of its moments (Bird & 
Beers, 1993). By definition, skewness, kurtosis and Al are 
equal to 0 for a gaussian dish, while TI to 1. In Table ^ 
we present the results; significance levels have been es¬ 
timated from Table 2 in Bird & Beers 1993. While the 
values obtained for skewness and Al cannot allow to re¬ 
ject the Gaussian hypothesis (significance level > 10%), 
both kurtosis and TI indicate departure from a Gaussian 


Fig. 6. Stripe density plot of radial velocities for the 
125 members of Abell 521. 



Fig. 7. Same as Fig. |^, but now the best fit Gaussians 
found by KMM for a three group partition is super¬ 
imposed and we have used a binning of 500 km/s. 


distribution at better than 10% significance level (Beers et 
al. 1991). This indicates that the dataset has more weight 
in the tails than a Gaussian of the same dispersion. 

As departure from normality and high values of veloc¬ 
ity dispersion can result from a mixing of several velocity 
distributions of smaller velocity dispersion with different 
locations, we have investigated various tests for the exis¬ 
tence of substructure in the cluster velocity distribution. 

We have therefore addressed the presence of gaps 
which can be a signature of sub-clustering (Beers et al. 
1991). Five significant gaps in the ordered velocity dataset 
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Table 1. ID substructure indicators for the 125 objects with 
quality flag=l in our dataset 


Indicator 

Value 

Significance 

AI 

0.470 

<0.20 

TI 

1.113 

<0.10 

Skewness 

0.087 

>0.20 

Kurtosis 

0.779 

<0.05 


Table 2. Weighted gaps in the data 


Velocity [km/s] 

Size 

Significance 

73455.4 

2.298 

0.030 

75243.4 

2.323 

0.030 

73704.0 

2.454 

0.030 

76738.1 

2.577 

0.014 

72150.0 

3.089 

0.002 


Table 3. 1-D statistical tests performed in ROSTAT package 
that exclude the hypothesis of a single gaussian distribution. In 
Cols. 1 and 2 we report the name and the value of the statistics, 
while Col. 3 indicate their significance levels. 


Statistical Test Value Significance 

a 

0.743 

<0.10 

w 

0.969 

0.08 

B2 

3.779 

0.04 

DIP 

0.023 

<0.05 

KS 

0.873 

0.10 

V 

1.707 

0.01 


0.162 

0.02 


0.162 

0.01 

A^ 

0.989 

0.01 



Fig. 8. Projected coordinates of the galaxies as¬ 
signed to the KMM partition. From top to bottom, 
KMM-A, KMM-B, KMM-C. The subgroups of the 
cluster identified by the isodensity contours in Figs. 
and ^ have been schematically represented (in next 
figures too). 


were detected. Fig. || shows the stripe density plots of ra¬ 
dial velocities of the 125 cluster galaxies and we have in¬ 
dicated the gap positions with an arrow, while in Table 
one can hnd the velocity of the object preceding the gap, 
the normalized size (i.e., the “importance”) of the gap it¬ 
self, and the probability of finding a normalized gap of this 
size and with the same position in a normal distribution. 
Two very significant gaps are detected respectively at ~ 
76740 km/s (probability lower than 1.4%) and ~ 72150 
km/s (probability lower than 0.2%). 

In addition, 9 of the 13 one-dimensional statistical 
tests of gaussianity performed by ROSTAT exclude the 
hypothesis of a single Gaussian distribution at better than 
10% significance level (see Table Among them, the Dip 
test is particularly indicative (Hartigan & Hartigan 1985). 
This tool tests the hypothesis that a sample is drawn from 
a unimodal, not necessarily Gaussian parent population. 
In the present case it rejects the unimodal hypothesis at a 
signihcance level better than 5%, confirming our previous 
results. 
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3.2. Partitioning the distribution in velocity space 

In order to separate possible velocity groups within our ve¬ 
locity dataset we have used the KMM mixture modeling 
algorithm of McLachlan & Basford (1988). This method 
has been shown to be very useful for detecting bimodal¬ 
ity in astronomical datasets (Ashman and Bird 1994), and 
can even be applied to detect multimodality. One of the 
major uncertainties however is the best choice of the num¬ 
ber of groups for the partition. Given the appearance of 
the velocity histogram and of the strip density plot, and 
the presence of two highly significant gaps, we have chosen 
as a first guess to fit three velocity groups around the mean 
velocities 70000, 74000 and 78500 km/s. The KMM algo¬ 
rithm then fits a 3-group partition from this guess and op¬ 
timizes the mean velocity for each group. However, in the 
case of a multimodal partition, the significance level is not 
determined accurately by the algorithm. The estimated P- 
value of 3 % suggests a strong rejection of the null hypoth¬ 
esis (unimodal gaussian distribution) but has to be taken 
only as a guideline. The low-velocity group A is best fitted 
by a Gaussian with parameters: Cbi = 71127tjgg km/s 
and Sbi = 678+gg km/s (17 objects). The major group 
B (103 objects) is found with Gbi = 74249/jgQ km/s and 
Sbi = 879t®i km/s. At the high velocity tail, 5 galaxies 
are found, with location Gbi = 78416/ji84 km/s, populat¬ 
ing group C. The three gaussians corresponding to these 
partitions are displayed in Fig. ^ 

We present in Fig. ^ the projected positions of the 
galaxies assigned to the three partitions. The position of 
the main overdensities emerged in the projected density 
maps are also displayed for an easier reading. Galaxies 
in KMM-B are following the cross-like general pattern of 
the cluster. Half of the galaxies in KMM-A are located in 
the central ridge. Galaxies belonging to KMM-G are all 
located in the South-East region of the cluster. We have 
also tried partitions with a higher number of groups (re¬ 
spectively 4 and 5) which also show strong rejection of the 
unimodal hypothesis. However, as suggested by Ashman 
and Bird 1994, we have followed the Occam’s razor and 
selected the partition with the smaller number of groups 
(tri-partition). 


4. Combined analysis of sub-clustering in 
velocity/2D space 

4.1. Tomography of Abell 521 

As indicated by the moments analysis, the tails of the 
velocity distribution appear to be widely populated, sug¬ 
gesting the presence of interlopers. Moreover, the analysis 
of the number density maps of the cluster (Arnaud et al. 
2000 and Fig. of the present paper) shows a complex 
structure, with six groups in projected coordinates (gl 
to g6) along the main SW/NE direction of the cluster, a 
group around the BGG and a high density ridge in the 
direction perpendicular to the main axis of the cluster. 


In Fig. the values of the radial velocities (bottom), 
the mean radial velocities (middle), and the velocity dis¬ 
persions (top) are plotted as a function of an angular ra¬ 
dius; velocity location and scale have been computed in 
concentric shells with a fixed number of objects (20), in 
order to have comparable statistics. Different colors in the 
velocity vs. radius plots have been used to visualize the 
different groups identified on the isodensity map. 

For the left panels, the center is taken at the position of 
the main X-ray cluster (Arnaud et al. 2000), which roughly 
coincides with the barycentre of the galaxy distribution. 
For the plots of the central column the origin is centered 
on the BGG position, while in the the right panels we 
used the projected coordinate along the cluster main axis 
S2 (North-West/South-East). Negative values correspond 
to the South-East extremity of the cluster, zero to the 
center, and positive values to the North-West extent. 

In Fig. ^ galaxies have been circled with different 
colors corresponding to different velocity bins. Fig. ^ and 
Fig. ^ can then be analyzed together to understand the 
variation of the velocity distribution in the field. The 
galaxies belonging to the region of the “ridge” SI, cen¬ 
tered on the barycentre position and extending perpendic¬ 
ularly to the main axis S2, are color-encoded as blue open 
squares in Fig. For this particular region, we obtain sys¬ 
tematically a lower mean velocity and a higher dispersion 
than for the whole cluster. This is particularly apparent in 
Fig. ^ top and medium, right panels. This region consists 
of several clumps in projected density map, but galaxies 
in the whole velocity range [70000-78000] km/s populate 
the various clumps. However, it is the large number of low 
velocity objects that lowers the value of location in this 
region as compared to that obtained for the whole cluster. 

At ^ 1.5 arcmin from the barycentre of the cluster, 
in the North-East direction, we find a compact region of 
galaxies at similar velocities corresponding to the BGG 
region (galaxies color-encoded as red squares in the bot¬ 
tom panels of Fig. |^), which shows a significantly lower 
velocity dispersion and a slightly higher value of the mean 
velocity (top and medium panels of Fig. ^). 

A high velocity group of four galaxies (~ 78500 km/s) 
is detected at large radius (~ 4 arcmin) from the barycen¬ 
tre of the cluster (bottom panel of Fig. ^), with a gap 
in velocity of more than 1500 km/s as compared to oth¬ 
ers galaxies at the same radius. Thus, these objects are 
likely to be unbound to the cluster; they are located in the 
South-East extremity of the cluster (bottom, right panel 
of Fig. H), and three of them in the region defined as g2 
(green open squares in Fig. ||). This latter is the result of 
the superimposition along the line of sight of these three 
high velocity galaxies and a concentration of objects in 
the velocity range ([72000-74000] km/s). 

Several other clumps have a more homogeneous veloc¬ 
ity composition: gl (cyan open squares in Fig. ^), g3 (pur¬ 
ple open squares in Fig. ^ and g4 ( yellow open squares 
in Fig. gl is mainly populated with galaxies in the 
[74000-76000] km/s velocity range, while the greatest part 
of objects in g3 have velocity in the [74000-75000] km/s 
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Fig. 9. Left (from bottom to top): radial velocities of the 125 galaxies of our spectroscopic catalogue, 
and velocity location and scale of concentric shells with a fixed number of objects (20), centered on the 
optical barycentre of the cluster. Center: as before, but the new origin of the x-axis is the BCG position. 
Right: as before, but Xproj is the projected coordinate along the main axis of the cluster (NW/SE) 
and the center position is again on the optical barycentre of the cluster. Continuous lines represent the 
whole sample velocity scale (top) and location (mean and bottom). Dashed lines (bottom) indicate the 
[v — (Tv, V -b (Tv] interval. Different colors correspond to the different subgroups detected on the isodensity 
map: red open squares: BCG group - blue open squares: ridge region - cyan open squares: clump 
gl - green open squares: g2 - purple open squares: g3 - yellow open squares: g4 - green full 
squares: g5 - cyan full squares: g6 - black open squares: remaining objects 


bin. In particular, as one can see comparing the first and 
fifth bins of the central, top panels of Fig. g3 has a 
mean velocity and velocity dispersion quite close to the 
value of the BCG group. The North-East g4 clump also 
shows a mean velocity location comparable to the BCG 
region, and a very small value of velocity dispersion. 

Unfortunately, the number of measured radial veloc¬ 
ities within the previous mentioned substructures is not 


always large enough to obtain meaningful dynamical in¬ 
formation for each of them; we have thus divided the clus¬ 
ter into three regions with a number of objects sufficient 
to derive stable estimators of location and scale and to get 
enough statistics in velocity histograms without degrading 
too much the binning. These regions have been defined as 
three slices perpendicular to the main direction of the clus¬ 
ter (see Fig. 10). This choice has been motivated by the 
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Fig. 10. In black, galaxy iso-density contours for B<27 and I<20; the lowest iso-density contour corresponds to Icr 
level above the mean density in the field, the contours are spaced by O.Scr. The division in 3 slices has been shown. 
Circles correspond to the galaxies in our high quality velocity dataset; different colors correspond to different ranges 
in radial velocity: blue: v^. = 69000-^72000 - green: v^ = 72000-F73000 - cyan: v^ = 73000-^74000 - yellow: 
V,. = 74000^75000 - red: v,. = 75000^76000 purple: v^ = 76000^80000 


following reasons: first of all we wish to test separately the 
velocity distribution within the high density ridge SI per¬ 
pendicular to the direction of the main cluster S2, which 
was shown in previous analysis to be of specific interest. 
We therefore design the central slice (2) to include it. We 
also want to investigate the possible differences in velocity 
distribution between the North-West and South-East re¬ 
gions suggested in Fig. ||. The southern slice (1) includes 
the groups gl, g2 and g3; the northern slices (3) includes 
both the BCG region as well as the northern extensions 
embedding g4, g5 and g6 subgroups. In Fig. ^the velocity 
distributions for each slice are shown; the corresponding 
values for velocity location and scale have been reported 
in Table 

In the southern region (slice I), the mean velocity of 
the main structure, obtained when excluding the four high 
velocity galaxies, is Cbi = 73886(^g5®® km/s, and its scale 
is Sbi = 1117^86° km/s. The bimodal appearance of 

the main structure histogram motivated us to try a KMM 
partitioning of the distribution. A very good agreement 
was found (0.929 significance level) for a fit by two gaus- 
sians centered respectively at 73070 km/s and 75140 km/s, 
and with velocity dispersions of Sbi = 5701*2 km/s and 
Sbi = 4981*^ km/s. This indicates that in the region 


there is probably a mix of two kinematically distinct pop¬ 
ulations in addition to the high velocity group, and reflects 
the previous mentioned difference in velocity distribution 
within the clumps gl and g3 with respect to g2. 

In slice 2, corresponding to the ridge region, the ve¬ 
locity distribution shows a very dispersed boxy shape, 
with a velocity dispersion reaching the very large value of 
1780;ii^2 km/s, in good agreement with previous values 
by Maurogordato et al. 2000. In agreement with previous 
results, the location of the central slice is the lowest one 
(Cbi = 736251*^^ km/s). 

The northern region (slice 3) shows a higher location 
than the previous slices, Cbi = 74300lj()4 km/s, which is 
comparable to the velocity of the BCG (74357±44 km/s), 
and a lower velocity dispersion (Sbi = 839^134 km/s). 
The shape of the distribution is symmetrical, although 
few galaxies are still present in the low velocity tail. We 
have addressed the dynamics of the region immediately 
surrounding the BCG; this galaxy has a complex struc¬ 
ture, with a series of bright knots embedded in a lower 
density arclike structure at 24 hys”^ kpc of its center 
(Maurogordato et al. 1996, Maurogordato et al. 2000). 
We have measured the redshifts of the four bright knots 
in the arclike structure (see Table ||), showing that these 
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objects belong to the cluster and are not gravitation¬ 
ally lensed background objects. The central location in a 
- 240 h75-i kpc region around the BCG (ten objects with 
measured velocities including the BCG and its multiple 
nuclei) is Cbi = 74340l^g2 km/s. This values is very close 
to the BCG radial velocity (74357±44 km/s); the velocity 
scale is very low, 256l®33 km/s, and of the same order 
as the BCG internal velocity dispersion, 368±463 km/s 
(Maurogordato et al. 2000) strongly suggesting it is bound 
to the BCG. 

When comparing the velocity distribution of the three 
slices, the northern region clearly shows a higher location 
and a lower velocity dispersion than the other regions, 
while a very high velocity dispersion is observed within 
the central ridge. These results confirm the trends that 
emerged in the velocity profiles in Fig. |^. 
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4.2. Kinematical indicators of sub-clustering 

Our previous analysis has clearly shown the presence of 
sub-clustering in the projected density distribution and 
the departure from gaussianity of the velocity distribution, 
indicating that the system has not yet reached equilib¬ 
rium. As a second step, we have performed a more system¬ 
atic search of sub-structures by addressing directly corre¬ 
lated deviations in position and velocity distributions. We 
have applied several classical methods that quantify the 
amount of substructures in galaxy clusters using positions 
and velocities. 

In Table || we list the actual values for A (Dressier & 
Shectman, 1988), e (Bird, 1994), and a (West & Bothun, 
1990) parameters and the significance of the corresponding 
tests, obtained through the bootstrap technique and by 
normalizing with 1000 Monte Garlo simulations. 

Assuming that these tests reject the null hypothesis 
if the significance level is less than 10%, only the A test 
finds evidences of subclustering at a high confidence level. 
This result is further investigated by using the kinematical 
estimators introduced by Girardi et al. (1997), which take 
into account separately the departures of the local mean 
(Sy) and dispersion (Js) from the global measurements for 
the whole cluster. Low values of local velocity dispersion 
will give high values of dsj while dy will be high in case of 
strong departures (i.e. more than one a) of the local mean 
velocity with respect to the global mean. In Fig. |l^ each 
galaxy is represented by a circle whose diameter is propor¬ 
tional to (top panel), (middle panel) and (bot¬ 
tom panel). The top panel shows several areas with many 
large circles, which indicate correlated spatial and kine¬ 
matic variations. In the middle panel, very large circles 
are present in the South-East region of the cluster, corre¬ 
sponding to the group of background galaxies detected in 
previous sections. One can also note several large circles 
in the ridge region, corresponding to the presence of the 
previously detected low velocity group. Moreover, in the 
bottom panel, we detect three clumps in which the veloc¬ 
ity dispersion is effectively low; one is in the South-East 
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Fig. 11. Velocity distributions of the three slices of 
Fig. 0; the contribution due to different morphologi¬ 
cal types has been shown using different colors (early: 
red, late: blue, undetermined color objects: black 
shading). A binning of 1000 km/s has been used. 
Velocity location and scale found with ROSTAT are 
also shown. 
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Fig. 12. Projected positions of the galaxies in our 
spectroscopic sample, represented by circles with di¬ 
mensions weighted by the estimators of Dressier & 
Shectman 1988 (A, top panel) and Girardi et al. 
1997 (Jy, middle panel , and ds, bottom panel). In 
the top panel, concentrations of large circles indicate 
a correlated spatial and kinematical variation. In the 
middle panel, strong deviations of the local velocity 
as compared to the global one, while in the bottom 
panel, a low value of the local velocity dispersion as 
compared to the mean one. Crosses correspond to 
circles larger than the box (see text for details). 


Table 4. Velocity distribution properties of the various sub¬ 
samples of A521 galaxies 


Subsample 

Galaxy nb. 

Cbi 

[km/s] 

Sbi 

[km/s] 

Whole sample 

125 

740191“" 

13251)“ 

KMM-A 

17 

71127^111 

6781?? 

KMM-B 

103 

742491^^ 

8791?) 

KMM-C 

5 

784161?“ 

- 

Non-emission line 
galaxies 

110 

74166l??2 

10871)“ 

Emission line 
galaxies 

15 

72390I923 

22501?“ 

Early 

72 

740831)“ 

11051)?® 

Late 

41 

734021“® 

17231??) 

Bright 

(Iab < Iab*) 

9 

743211)“ 

4 QI+265 

Intermediate 
(Iab* < Iab < lAB*-k2) 

76 

740001)“ 

13311)?? 

Faint 

(Iab > lAB*-k2) 

28 

735151)?? 

16351“) 

Slice 1 

44 

740371“? 

14541??? 

Slice 1 without 
the background group 

40 

738861)?? 

11171)?° 

Slice 2 

30 

736251?“ 

17801??) 

Slice 3 

51 

742981))? 

8391?)? 


Table 5. Radial velocities of the four blobs embedded in the 
arclike structure at 24 hrs”^ kpc of the BCG center. The names 
of the blobs refer to Maurogordato et al. 2000, Fig. 8. 


Blob 

Radial Velocity (km/s) 

Error on Velocity (km/s) 

B 

74340 

100 

C 

74341 

80 

D 

74325 

52 

E 

74205 

106 


region and corresponds to the structure identified as g3 
in our iso-density maps, while the second is centered on 
the BCG complex. The third clump consists of a group 
of four galaxies in the g4 region; it has such a low local 
velocity dispersion that the corresponding values are 
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Table 6. 3-D substructure indicators for the sample of 125 Table 7. Properties of the significant groups from the htree 
objects with quality flag=l in our dataset analysis 


Indicator 

Value 

Significance 

A 

164.8 

0.080 

e 

1.91x10+^'^ kg 

0.217 

a 

0.183 Mpc 

0.381 


larger than the box size of Fig. |l^ Therefore, we have 
represented these objects with circles of diameter equal 
to their (5s parameters, and we have indicated their pro¬ 
jected positions with crosses. Their radial velocities are in 
the range [74200-74500] km/s, as can clearly be seen in 
Fig. where the four galaxies are represented by yellow 
squares. We also note that in the region of the “ridge” 
and in the extreme South-East region, the dimension of 
the circles is the smallest, indicating in these regions the 
highest values of the velocity dispersion are reached. 

We have also applied the algorithm developed by Serna 
& Gerbal (1996) which identifies subclusters on the basis 
of dynamical arguments. It uses the hierarchical clustering 
algorithm to associate galaxies according to their relative 
binding energies. Fig. shows the map in projected coor¬ 
dinates visualizing the groups resulting from the substruc¬ 
ture analysis at the various levels. At the first level, the al¬ 
gorithm separates the data into the main cluster G1 of 120 
objects (open circles -I- stars) with Vmean = 73835 km/s 
and a = 1204 km/s, and a high velocity set of 5 ob¬ 
jects G2 (large filled points) in the southern region, with 
Vmean = 78418 km/s and cr = 500 km/s, previously identi¬ 
fied as KMM-C. The algorithm then splits again the main 
cluster G1 in two components of very different mass ra¬ 
tios: the main one Gil (93 objects, open circles) with a 
mean velocity similar to the value obtained for G1 , but 
now a velocity dispersion greatly reduced (tr = 940 km/s), 
and a low mass group G12 of five objects at higher ve¬ 
locity (75730 km/s) plotted as stars in the southern re¬ 
gion. These objects were already detected in the analy¬ 
sis of slice 1, as a possible higher velocity population of 
the southern region. At the last level of structure iden¬ 
tification, the algorithm finds out two significant groups: 
the first one at North (Gill), corresponding to the re¬ 
gion surrounding the BGG galaxy (squares), at a slightly 
higher mean velocity than the mean component (74290 
km/s) and presenting a low value of the velocity disper¬ 
sion (442 km/s), and a compact group southern of the 
ridge (circles), corresponding to g3 in previous Sects. The 
energy levels of the various groups are also provided (not 
displayed): the deepest energy levels at the bottom of the 
energy well corresponds to the BCG complex. The lower 
level is occupied by the BCG itself associated to blob A, 
joining with the three other blobs C, D and E (see Table H) 
at a slightly higher level. Associating to the BCG com¬ 
plex various galaxy pairs at low energy levels results into 
a bound central group around the BCG galaxy. It is in¬ 
teresting to note that the BGG group and g3 are detected 


group 

^mean 

[km/s] 

cr 

[km/s] 

(10^®M©) 

Ngai 

main 

74018 

1386 

1.96 

125 

G1 

73835 

1204 

1.64 

120 

G2 

78418 

504 

0.2 

5 

Gll 

73965 

930 

1.1 

93 

G12 

75730 

121 

0.006 

5 

Gill 

74290 

442 

0.045 

12 

G112 

74068 

574 

0.063 

7 


both by the dynamical estimator of Girardi et al. 1997, 
and by the htree algorithm. 

4.3. Dynamics of the groups 

As seen before, several groups are revealed from the sub¬ 
structure analysis. We have then performed the gravita¬ 
tional bound check as in the two-body problem (Beers, 
Geller and Huchra 1982): 

Vr^Rp < 2 G M sin^o; cosa, (1) 

where Vr is the relative line-of-sight velocity of the two 
considered clumps, Rp is the projected separation of the 
clump centers, a is the projected angle measured from the 
sky plane, and M is the total mass of the system. 

In a first step, we have applied it to the groups G1 and 
G2 previously defined by the h-tree method. The result 
is displayed in Fig. 0 , which shows that the system is 
unbound for nearly all values of the projection angle a. 
We have then eliminated the group G2 from the analysis 
as a background system of galaxies, and tested if the group 
Gill (corresponding to the system surrounding the BGG) 
and the group G112 (corresponding to g3) are bound to 
the remaining cluster. We considered each time a two- 
body problem, with one system being the tested group, 
and the other the remaining cluster. It results that the 
groups Gill and G112 are bound to the cluster for nearly 
all values of the projection angle a (4 to 90 degrees). 

5. Variation of dynamical properties with color 
and luminosity 

In this Sect, we use the color information to define early 
and late type galaxy subsamples, using the result that 
the bulk of early-type galaxies in all rich clusters usu¬ 
ally lie along a linear color-magnitude relation. This so 
called “red sequence” has been interpreted as a clear in¬ 
dication that all rich clusters contain a core population of 
passively, evolving elliptical galaxies, coeval and formed 
at high redshift (Ellis et al. 1997, Kodama et al. 1999, 
Gladders et al. 1998). In Fig. ^the (B-I) versus I color- 
magnitude diagram (CMD) is shown for the objects in a 
field of 15x15 arcmin^ (^3.25x3.25 hrs”^ Mpc^) centered 
on the optical barycentre of the cluster. 
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Fig. 13. Substructure analysis by the h-tree method: the map in projected coordinates is shown, with different symbols 
corresponding to the different groups identified. 


Due to the fact that we have no magnitude information 
for 12 of the 125 objects of A521, our velocity/magnitude 
catalogue is made up of 113 galaxies. 

The red sequence of the cluster has been characterized 
(slope, intercept and width) by considering galaxies in a 
smaller field (6.8x6.8 arcmin^, ~1.5x1.5 hys^^Mpc^) in 
order to reduce the higher contamination in periphery by 
field objects. We have used an algorithm of linear regres¬ 
sion plus an iterated Scr clip, following and refining the 
method of Gladders et al. (1998) (detailed in Ferrari et 
al. in prep.). The red sequence can then be described by 
the linear equation (B — 1)^^= —0.033 Iab + 3.273 with 
a width of 0.183 (see Fig. |^). Galaxies with measured 
redshift that lie within the identified red sequence and the 
three reddest objects of the cluster in the color-magnitude 
diagram have been classified as early type objects, while 
the galaxies with bluer colors on this diagram have been 
assumed to correspond to later types. 

We have analyzed the galaxy velocity distribution of 
our sample as a function of this approximate early/late 
classification; histograms of velocity and distribution in 
projected coordinates of these objects are displayed in 
Fig. Velocity locations and scales of the two subsam¬ 


ples calculated using the biweight technique are shown in 
Table |. 

We find significantly different velocity dispersions for 
the two samples, with an extremely high value for the late 
type sub-sample, while quite close to the value of the main 
cluster component (KMM-B) for the early type one. 

A trend in location is also detected, with the late type 
objects located at a lower velocity as compared to the early 
type sample. A Kolmogorov-Smirnov test rejects the hy¬ 
pothesis that the two velocity datasets are drawn from the 
same parent population with a significance level of ^5.5% 
(see table ||). The distribution in projected coordinates 
is also very different for the two samples (Fig. 0 right 
panels). Most of the early type galaxies are aligned with 
the S2 axis and populate the region of the BCG and the 
groups, with some sparse objects in the outskirts, whereas 
the galaxies of the late type sample are mainly located 
along the SI filament and in the external regions of the 
cluster, with very few objects in the South. This is con¬ 
firmed by the velocity histograms plotted in Fig. for 
each slice of Fig. ^ where the counts in velocity bins are 
encoded with a different color for each morphological type 
(blue for late, red for early). The region of the ridge (slice 



14 


Ferrari, C. et al.: Multiple merging events in Abell 521 



Fig. 14. Boundary between the bound and unbound 
orbits for the main cluster G1 and the high veloc¬ 
ity system G2. The domain left of the curve corre¬ 
sponds to bound orbits, while the right part to un¬ 
bound ones. The vertical line corresponds to the ve¬ 
locity difference between the mean location of the 
two groups. The corresponding 68% confidence re¬ 
gion is shown with the cross-hatching. 



16 18 20 22 


Fig. 15. The color-magnitude diagram (B-I) vs. I for 
all the galaxies in a 15x15 arcmin^ area containing 
the core of the cluster. The objects represented with 
a circle have spectroscopic redshift and belong to the 
cluster. The solid straight line is the linear best-fit 
found for the red sequence of the elliptical galaxies 
of the cluster, while the dashed line represents the 
“boundary” between early and late type objects on 
the GMD. 


Fig. 16. Top: Velocity distribution with a binning of 
800 km/s (left) and projected coordinates (right) of 
early type galaxies. Bottom: as before, but for the 
late type subsample. 

2) is the richest in late type objects (~ 50% versus ^ 28% 
for slices 1 and 3), in particular in the low velocity tail. 

We have then investigated if this strong morphological 
segregation could be related at least partially to a lumi¬ 
nosity segregation. Using the previous classification, we 
have converted apparent magnitudes to absolute ones, as¬ 
suming k and e-corrections for E and Sa models given by 
Poggianti (1997). We have then divided the catalogue in 
three subsamples in absolute magnitudes, corresponding 
to Iab <Iab*, Iab* < Iab < Iab * + 2, Iab >Iab * + 2 (with 
Iab* = —22.7, derived using the i'* value in Goto et al. 
2002 and following the indication of Fukugita et al. 1995 
for transforming i’ in Iab), in order to test if the spa¬ 
tial and velocity distribution differs for various luminosity 
classes within the cluster. In table ^ we report the results 
for each sub-sample. 

In Fig. Ill, we display the projected distributions and 
the velocity histograms for the three tested subsamples. 
As a first evidence, the velocity dispersion increases dras¬ 
tically from brightest to faintest objects, with more than 
a factor two between the value obtained for the Iab<Iab* 
sample towards the Iab>Iab* + 2 one. There is also some 
indication of shift in location, as the brightest sub-sample 
is characterized by a mean velocity ~700 km/s higher than 
the fainter ones. 

The brightest galaxies lie only on the main axis of the 
cluster (S2), while the faintest ones are preferentially lo¬ 
cated along SI with some galaxies in the outskirts. We 
have then tested for correlations between the pseudo- 
spectral type and luminosity class distributions. The late 
type subsample and the faint one have extremely similar 
distributions both in velocity and in projected positions; 
in both cases, a K-S test indicates that they are drawn 
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Fig. 17. Spatial and velocity (with a binning of 800 
km/s) distribution for the three Iab magnitude sub¬ 
samples of A521 galaxies. Top: Iab<-22.7 - Middle: 
-22.7< Iab <-20.7 - Bottom: Iab>-20.7 


Table 8. Significances (%) given by the Kolmogorov-Smirnov 
Test for the velocity distributions of the three slices of Fig. 



Slice 1 

Slice 2 

Slice 3 

Slice 1 

- 

37.4 

27.1 

Slice 2 

- 

- 

12.3 


from the same distribution function with high significance 
levels (more than 99% and 94%, Table ||). On the other 
side, both the velocity and the projected position distribu¬ 
tions of the early type and of the intermediately luminous 
objects are drawn from the same parent population with 
a significance level of more than 97%, while the early type 
and the brightest galaxies have significance levels higher 
than 80%. 

We have also analyzed separately the population of 
emission lines (15 objects) versus non emission lines galax¬ 


Table 9. Significances given by the Kolmogorov-Smirnov Test 
for the velocity (top) and spatial (bottom) distributions of the 
various detected sub-samples 


Significances (%) 
for velocity distributions 


Early 

Late 

Bright 

Interm. 

Faint 

Early 

- 

5.7 

81.8 

97.9 

22.2 

Late 

- 

- 

34.8 

30.7 

99.7 

Bright 

- 

- 

- 

58.7 

20.4 

Interm. 

- 

- 

- 

- 

56.1 

Significances (%) 

for projected position distributions 


Early 

Late 

Bright 

Interm. 

Faint 

Early 

- 

47.4 

87.2 

99.9 

69.8 

Late 

- 

- 

77.0 

92.6 

94.2 

Bright 

- 

- 

- 

77.0 

64.3 

Interm. 

- 

- 

- 

- 

35.6 


ies (110 objects) (Fig. ^8|). When excluding emission lines 
galaxies, the location remains comparable but the velocity 
dispersion is notably smaller than that of the whole sam¬ 
ple (Table As shown by the histogram in Fig. the 
velocity distribution of emission lines galaxies consists of a 
major concentration at low velocity (around 70000 km/s) 
and several isolated objects spanning the whole range of 
the cluster. Considering these objects together with the 
non-emission lines objects results in strongly enhancing 
the global velocity dispersion. 

6. Discussion and Conclusions 

The velocity distribution of Abell 521 is definitively 
very complex. The large value of the velocity dispersion 
(^ 1325 km/s) of our whole spectroscopical sample of 125 
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Fig. 18. Left: Velocity distribution (with a binning of 500 km/s) for galaxies with (white) and without 
(shading) emission lines. Right: Projected coordinates of the emission line galaxies. Solid squares: 
Vr = 69000-F72000 - Open squares: Vr = 72000-^73000 - Solid circles: = 74000-^75000 - Stars: 
Yr = 75000-b76000 - Triangles: = 76000-^80000. On the electronic version of the article, colors of 

the symbols have the same meaning as in Fig. 10, 


galaxies clearly results from the mixture of several compo¬ 
nents. The main features of the velocity distribution are 
summarized hereafter: 

i) a high velocity tail in the velocity distribution is 
revealed by a KMM partition. From the analysis of the 
projected positions and of the velocity/radius diagram, 
these objects are found to lie in the South-East region 
at about 870 h 75 ~^ kpc from the X-ray main center of 
the cluster, and to have velocities higher than 1500 km/s 
as compared to the other objects at the same radius. The 
two-body criteria shows that the probability that this sys¬ 
tem is not bound to the main component of the cluster is 
quite high. These objects are probably field galaxies or a 
loose background group. When excluding these galaxies, 
the mean location remains unchanged (~ 74000 km/s), 
but the scale reduces slightly (^ 1200 km/s). 

ii) The velocity distribution in the high density central 
ridge shows a velocity location systematically lower than 
for the whole cluster, and a very high velocity dispersion 
(1780t?34 

km/s). Moreover the ratio of late/early type ob¬ 
jects are higher than in the other slices (0.95 versus 0.46 
in slice 1, and 0.48 in slice 3). This region is also particu¬ 
larly rich in emission-line galaxies, as it contains one third 
of the emission line objects detected in the whole sample. 
The late type, emission-line objects coincide with the low 
velocity tail of the distribution (v<72000 km/s. Figs, p] 
and [l^) . Therefore, we are witnessing a very unusual con¬ 
figuration in the core of the cluster: a filamentary structure 
of ~ 1 h 75 ~^ Mpc, with galaxies showing a very broad ve¬ 
locity dispersion with an excess of low velocity objects, 
colors typically bluer than the mean, and the presence of 
several emission line galaxies. 


hi) Various groups lie along the main NW/SE axis of 
the cluster. First, the region including the complex around 
the BCG (~ 240 h 75 “^ kpc) appears as a strongly bound 
system with a very low velocity dispersion (^ 250 km/s) 
typical of a group and a location higher than the whole 
cluster(~ 74340 km/s). Other groups with comparable 
values of location are observed: g4 at the NE extent and 
the southern group g3. The other groups are less well- 
sampled, so the following results have to be taken with 
caution. The northern group g6 shows a slightly lower lo¬ 
cation, while at South, gl seems to belong to a higher 
velocity complex, and g2 result of superposition effects. 

From these results, we can refine the scenario of forma¬ 
tion of the cluster. The northern region is characterized by 
a lower velocity dispersion and a slightly higher location; it 
hosts a group dynamically bound to the BCG; it is clearly 
associated to the compact group in X-ray, which is prob¬ 
ably falling on the main cluster (Arnaud et al. 2000). The 
small difference in the mean velocity of the northern re¬ 
gion as compared to the whole cluster (^ 250 km/s in the 
rest-frame of the cluster) suggests that the merging occurs 
partly in the plane of the sky, along the North-West/South 
East direction S2. This direction emerges as the main axis 
of the ongoing merging event. Most of the detected clumps 
are aligned along this direction. Moreover, the early type 
objects as well as the brightest ones (L> L*) follow the 
general NW/SE skeleton of the cluster. 

However, our scenario has also to reproduce the pe¬ 
culiar features that appear in the central region of the 
cluster, in particular its high velocity dispersion, its lower 
location, and its filamentary NE/SW structure. The high 
velocity dispersion is mostly due to the presence of the 
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low velocity component detected by both the KMM par¬ 
tition and the velocity profiles. It could consist of a fore¬ 
ground group currently interacting with the main cluster. 
It is separated in radial velocity space by ^ 3000 km/s 
from the main cluster component. This high velocity bulk 
flow could result from a recent merger with a significant 
component along the line of sight direction. Such high¬ 
speed encounters can be detected in merging clusters, as 
shown for instance in the case of C10024-I-1654 (Czoske et 
al. 2002). The high density ridge SI would then be the 
projection of the merging axis along the plane of the sky. 
The large fraction of late type and star-forming objects 
along SI, coincident with the objects of the low velocity 
partition KMM-A, corroborates this hypothesis. In fact, a 
strong compression of the gas perpendicularly to the axis 
of the merger between the two components is expected 
during the merging event, which could trigger star forma¬ 
tion in this region (Caldwell et al. 1993, Caldwell & Rose 
1997, Bekki 1999). Moreover, the brightest galaxy in the 
eastern side of the ridge shows the same orientation as SI, 
and a velocity of ~ 72000 km/s. This object could in fact 
be the original brightest galaxy of the group which has 
collided the main cluster. 

These results imply that Abell 521 is the outcome of 
multiple merger processes at various stages. A denser sam¬ 
pling in velocity of the various groups, combined with a 
wider angular coverage are planned in order to understand 
the large-scale dynamics of this particular cluster. 
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Table 10. Velocity data 


# 

RUN 

R.A. 

(2000) 

DEC. 

(2000) 

HEL.VEL. 

v(km/s) 

ERROR Quality flag Emission lines 

Av(km/s) 

1 

ES02 

04:53:47.91 

-10:11:23.52 

89027 

89 

3 

NO 

2 

ES02 

04:53:48.34 

-10:11:49.60 

73101 

39 

1 

NO 

3 

ESOl 

04:53:49.70 

-10:13:32.42 

33000 

187 

3 

NO 

4 

ES02 

04:53:49.94 

-10:11:18.81 

75683 

51 

1 

NO 

5 

ESOl 

04:53:50.26 

-10:12:19.00 

25724 

109 

1 

Ha,H/3,OIII 

6 

ES02 

04:53:50.96 

-10:10:17.04 

74317 

100 

3 

NO 

7 

ESOl 

04:53:51.55 

-10:12:42.90 

-2 

-2 

1 

NO 

8 

ES02 

04:53:52.01 

-10:11:52.27 

-2 

-2 

1 

NO 

9 

ESOl 

04:53:52.05 

-10:12:31.30 

74010 

95 

1 

NO 

10 

ESOl 

04:53:52.63 

-10:13:23.50 

-2 

-2 

1 

NO 

11 

ES02 

04:53:53.19 

-10:10:44.00 

-1 

-1 

4 

NO 

12 

ESOl 

04:53:53.86 

-10:12:23.32 

71668 

124 

2 

NO 

13 

ES02 

04:53:54.46 

-10:12:29.41 

-2 

-2 

1 

NO 

14 

ESOl 

04:53:54.75 

-10:12:19.65 

74144 

66 

1 

NO 

15 

ES02 

04:53:54.81 

-10:10:51.80 

-2 

-2 

1 

NO 

16 

ESOl 

04:53:55.45 

-10:15:01.01 

74017 

73 

1 

NO 

17 

ES02 

04:53:56.40 

-10:10:02.15 

74784 

57 

1 

NO 

18 

ESOl 

04:53:56.44 

-10:13:26.80 

76181 

139 

1 

NO 

19 

ESOl 

04:53:56.56 

-10:14:46.29 

75202 

122 

1 

NO 

20 

ES02 

04:53:56.66 

-10:11:41.18 

-1 

-1 

4 

NO 

21 

ESOl 

04:53:57.33 

-10:14:09.88 

67749 

86 

2 

NO 

22 

ESOl 

04:53:57.82 

-10:14:24.45 

-1 

-1 

4 

NO 

23 

ES02 

04:53:58.10 

-10:12:09.25 

74772 

63 

1 

NO 

24 

ESOl 

04:53:58.60 

-10:14:20.41 

74584 

80 

1 

NO 

25 

ESOl 

04:53:58.77 

-10:13:25.85 

71170 

71 

1 

NO 

26 

ES02 

04:53:58.95 

-10:11:26.08 

72694 

38 

1 

NO 

27 

ES02 

04:53:59.29 

-10:13:01.14 

121075 

78 

2 

NO 

28 

ESOl 

04:53:59.39 

-10:13:13.00 

-2 

-2 

1 

NO 

29 

ES02 

04:53:59.51 

-10:12:59.92 

73455 

34 

1 

NO 

30 

ESOl 

04:53:59.93 

-10:13:01.63 

73033 

52 

1 

NO 

31 

ES02 

04:54:00.17 

-10:11:14.58 

74250 

40 

1 

NO 

32 

ESOl 

04:54:00.27 

-10:14:05.13 

74574 

112 

1 

NO 

33 

ES02 

04:54:00.54 

-10:11:55.59 

76219 

49 

1 

NO 

34 

ESOl 

04:54:00.70 

-10:13:09.70 

182923 

129 

1 

Oil 

35 

ES02 

04:54:01.17 

-10:12:23.14 

74657 

114 

1 

NO 

36 

ES02 

04:54:01.41 

-10:11:16.30 

71490 

100 

1 

OII,OIIIa,b,H/3 

37 

ESOl 

04:54:01.62 

-10:14:07.20 

73928 

67 

1 

NO 

38 

ESOl 

04:54:01.75 

-10:12:34.47 

99494 

90 

1 

OII,H/? 

39 

ES02 

04:54:02.16 

-10:11:09.73 

193062 

100 

3 

NO 

40 

ES02 

04:54:02.26 

-10:12:49.86 

99810 

100 

1 

OII,Balmer,H/? 

41 

ESOl 

04:54:02.34 

-10:14:04.40 

70301 

101 

1 

NO 

42 

ESOl 

04:54:02.61 

-10:13:01.78 

73195 

92 

2 

Oil 

43 

ES02 

04:54:02.97 

-10:11:29.06 

-2 

-2 

1 

NO 

44 

ESOl 

04:54:03.26 

-10:13:22.40 

39765 

88 

1 

NO 

45 

ES02 

04:54:03.91 

-10:12:15.68 

74900 

111 

1 

NO 

46 

ESOl 

04:54:04.02 

-10:14:42.44 

73981 

98 

1 

NO 

47 

ES02 

04:54:04.02 

-10:09:22.45 

-2 

-2 

1 

NO 

48 

ESOl 

04:54:04.12 

-10:14:21.10 

-2 

-2 

1 

NO 

49 

ES02 

04:54:04.48 

-10:11:54.81 

-2 

-2 

1 

NO 

50 

ESOl 

04:54:04.58 

-10:14:14.43 

218349 

70 

2 

NO 
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# 

RUN 

R.A. 

(2000) 

DEC. 

(2000) 

HEL.VEL 

v(km/s) 

ERROR 

Av(km/ 

Quality flag Emission lines 

s) 

51 

ES02 

04:54:04.63 

-10:17:31.71 

72713 

62 

1 

NO 

52 

ESOl 

04:54:04.90 

-10:12:03.68 

72802 

82 

1 

OII,H/3,OIIIa,b 

53 

ESOl 

04:54:05.06 

-10:16:04.27 

69854 

89 

1 

H/?,OIIIa,b 

54 

ES02 

04:54:05.24 

-10:10:31.13 

1 

0 

1 

QSO 

55 

ESOl 

04:54:05.48 

-10:14:10.80 

74959 

112 

1 

NO 

56 

ES02 

04:54:05.55 

-10:17:37.33 

74862 

71 

1 

NO 

57 

ESOl 

04:54:05.60 

-10:11:15.50 

-2 

-2 

1 

NO 

58 

ES02 

04:54:05.79 

-10:11:41.71 

74225 

77 

1 

NO 

59 

ESOl 

04:54:05.98 

-10:13:16.60 

74571 

50 

1 

NO 

60 

ES02 

04:54:06.20 

-10:12:39.05 

-2 

-2 

1 

NO 

61 

ES02 

04:54:06.36 

-10:10:49.94 

74214 

83 

1 

OII,OIIIb 

62 

ES02 

04:54:06.36 

-10:18:14.36 

76171 

59 

1 

NO 

63 

ESOl 

04:54:06.37 

-10:15:59.98 

109406 

194 

1 

OII,H/3 

64 

ESOl 

04:54:06.43 

-10:13:24.76 

74340 

100 

2 

NO 

65 

ESOl 

04:54:06.56 

-10:13:21.32 

74341 

80 

2 

NO 

66 

ESOl 

04:54:06.64 

-10:15:36.00 

72394 

66 

2 

NO 

67 

ES02 

04:54:06.65 

-10:11:06.49 

-1 

-1 

4 

NO 

68 

ESOl 

04:54:06.94 

-10:11:46.43 

-2 

-2 

1 

NO 

69 

ESOl 

04:54:07.02 

-10:16:04.20 

44219 

121 

2 

NO 

70 

ESOl 

04:54:07.03 

-10:12:07.51 

74475 

79 

1 

NO 

71 

ESOl 

04:54:07.06 

-10:17:55.90 

74343 

107 

1 

NO 

72 

ES02 

04:54:07.10 

-10:13:16.76 

74205 

106 

1 

NO 

73 

ESOl 

04:54:07.12 

-10:12:41.25 

110097 

121 

1 

OII,H/3,OIIIa,b 

74 

ESOl 

04:54:07.14 

-10:15:11.36 

74110 

77 

1 

NO 

75 

ES02 

04:54:07.97 

-10:17:41.83 

-1 

-1 

4 

NO 

76 

ESOl 

04:54:08.05 

-10:14:01.90 

74334 

70 

1 

NO 

77 

ES02 

04:54:08.11 

-10:10:20.47 

74329 

43 

1 

NO 

78 

ES02 

04:54:08.14 

-10:11:11.79 

74536 

79 

1 

NO 

79 

ESOl 

04:54:08.23 

-10:14:58.59 

68502 

69 

3 

NO 

80 

ES02 

04:54:08.28 

-10:14:32.30 

74369 

100 

2 

NO 

81 

ESOl 

04:54:08.31 

-10:17:13.59 

55952 

37 

2 

H/? 

82 

ESOl 

04:54:08.41 

-10:12:41.95 

75867 

77 

1 

NO 

83 

ESOl 

04:54:08.56 

-10:15:13.65 

102858 

116 

3 

NO 

84 

ESOl 

04:54:08.60 

-10:15:50.10 

75243 

62 

1 

NO 

85 

ES02 

04:54:08.71 

-10:18:19.19 

71848 

42 

1 

NO 

86 

ES02 

04:54:08.83 

-10:10:48.42 

74259 

48 

1 

NO 

87 

ESOl 

04:54:08.88 

-10:14:58.00 

73236 

64 

1 

NO 

88 

ES02 

04:54:08.99 

-10:11:06.96 

109914 

100 

1 

on 

89 

ESOl 

04:54:09.08 

-10:15:34.90 

74438 

81 

1 

NO 

90 

ESOl 

04:54:09.33 

-10:15:10.08 

40747 

66 

3 

NO 

91 

ES02 

04:54:09.37 

-10:14:10.47 

76738 

85 

1 

NO 

92 

ESOl 

04:54:09.47 

-10:17:11.80 

75526 

51 

1 

NO 

93 

ES02 

04:54:09.50 

-10:16:44.87 

76202 

114 

2 

NO 

94 

ESOl 

04:54:09.86 

-10:16:03.13 

72919 

102 

1 

NO 

95 

ESOl 

04:54:09.90 

-10:14:13.40 

72972 

47 

1 

NO 

96 

ESOl 

04:54:09.90 

-10:17:34.20 

78810 

89 

1 

NO 

97 

ESOl 

04:54:10.15 

-10:12:55.40 

-2 

-2 

1 

NO 

98 

ESOl 

04:54:10.39 

-10:14:46.00 

87413 

90 

1 

NO 

99 

ESOl 

04:54:10.58 

-10:15:10.54 

106523 

114 

1 

OIIIa,b,Ha 

100 

ES02 

04:54:10.72 

-10:11:50.05 

-1 

-1 

4 

NO 
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# 

RUN 

R.A. 

(2000) 

DEC. 

(2000) 

HEL.VEL 

v(km/s) 

.. ERROR 
Av(km/ 

Quality flag Emission lines 

s) 

101 

ES02 

04:54:11.06 

-10:17:35.16 

75862 

50 

1 

NO 

102 

ESOl 

04:54:11.07 

-10:15:32.80 

-2 

-2 

1 

NO 

103 

ES02 

04:54:11.54 

-10:18:12.23 

-2 

-2 

1 

NO 

104 

ESOl 

04:54:11.57 

-10:15:13.02 

35443 

50 

2 

NO 

105 

ES02 

04:54:11.72 

-10:14:35.28 

70952 

100 

1 

OII,OIIIa,b,H/? 

106 

ES02 

04:54:12.02 

-10:10:27.20 

-2 

-2 

1 

NO 

107 

ESOl 

04:54:12.02 

-10:17:59.30 

-2 

-2 

1 

NO 

108 

ES02 

04:54:12.31 

-10:16:33.51 

-2 

-2 

1 

NO 

109 

ES02 

04:54:12.39 

-10:14:13.68 

72076 

60 

1 

NO 

no 

ES02 

04:54:12.50 

-10:14:20.31 

70102 

15 

1 

OII,OIIIa,b,H/? 

111 

ESOl 

04:54:12.73 

-10:15:51.20 

73965 

36 

1 

NO 

112 

ESOl 

04:54:13.02 

-10:15:51.80 

72985 

143 

1 

NO 

113 

ESOl 

04:54:13.09 

-10:13:52.61 

72067 

100 

1 

NO 

114 

ES02 

04:54:13.10 

-10:14:17.53 

36789 

100 

3 

NO 

115 

ESOl 

04:54:13.14 

-10:13:21.39 

28080 

59 

1 

NO 

116 

ESOl 

04:54:13.19 

-10:16:08.94 

93150 

120 

2 

NO 

117 

ESOl 

04:54:13.37 

-10:13:58.40 

94388 

113 

2 

NO 

118 

ESOl 

04:54:13.39 

-10:15:17.74 

-1 

-1 

4 

NO 

119 

ES02 

04:54:13.44 

-10:11:17.66 

-2 

-2 

1 

NO 

120 

ES02 

04:54:13.54 

-10:17:27.27 

131100 

100 

2 

on 

121 

ES02 

04:54:13.90 

-10:13:32.32 

88714 

100 

1 

OIIjBalmer 

122 

ES02 

04:54:13.97 

-10:11:13.55 

124290 

100 

3 

NO 

123 

ESOl 

04:54:13.97 

-10:18:02.64 

76000 

99 

2 

NO 

124 

ES02 

04:54:14.10 

-10:14:07.49 

106911 

89 

2 

NO 

125 

ESOl 

04:54:14.21 

-10:13:43.40 

107079 

98 

1 

NO 

126 

ESOl 

04:54:14.26 

-10:13:54.26 

88284 

94 

2 

NO 

127 

ES02 

04:54:14.53 

-10:14:43.21 

89310 

100 

2 

OII,H/3,OIIIb 

128 

ESOl 

04:54:14.55 

-10:12:35.33 

73619 

144 

2 

NO 

129 

ESOl 

04:54:14.69 

-10:18:23.59 

57498 

70 

3 

NO 

130 

ES02 

04:54:14.77 

-10:17:47.29 

70191 

100 

1 

OII,OIIIa,b,H/? 

131 

ES02 

04:54:14.99 

-10:11:21.55 

86628 

54 

1 

on 

132 

ESOl 

04:54:15.03 

-10:13:48.87 

127760 

131 

1 

NO 

133 

ESOl 

04:54:15.48 

-10:13:54.17 

75229 

90 

1 

NO 

134 

ES02 

04:54:15.62 

-10:17:20.94 

130800 

100 

1 

OIIjBalmer 

135 

ES02 

04:54:15.68 

-10:10:09.49 

109765 

40 

1 

NO 

136 

ESOl 

04:54:15.76 

-10:14:01.20 

49153 

124 

2 

NO 

137 

ESOl 

04:54:15.81 

-10:16:47.36 

71874 

101 

1 

NO 

138 

ESOl 

04:54:15.88 

-10:15:05.79 

78411 

55 

2 

NO 

139 

ESOl 

04:54:16.01 

-10:16:11.60 

74932 

92 

1 

NO 

140 

ESOl 

04:54:16.05 

-10:12:57.66 

72704 

93 

1 

NO 

141 

ES02 

04:54:16.08 

-10:10:08.86 

73830 

122 

1 

NO 

142 

ES02 

04:54:16.13 

-10:12:08.33 

-2 

-2 

1 

NO 

143 

ESOl 

04:54:16.34 

-10:16:04.61 

74282 

61 

1 

NO 

144 

ESOl 

04:54:16.44 

-10:15:09.60 

73398 

no 

1 

NO 

145 

ES02 

04:54:16.52 

-10:17:54.33 

107203 

84 

2 

NO 

146 

ESOl 

04:54:16.60 

-10:14:30.60 

73428 

103 

2 

NO 

147 

ES02 

04:54:16.92 

-10:18:09.91 

73070 

31 

1 

NO 

148 

ESOl 

04:54:17.12 

-10:14:25.17 

71017 

130 

2 

NO 

149 

ES02 

04:54:17.38 

-10:10:57.66 

74409 

76 

1 

NO 

150 

ES02 

04:54:17.54 

-10:13:05.01 

98004 

108 

2 

NO 



Ferrari, C. et al.: Multiple merging events in Abell 521 


21 


# 

RUN 

R.A. 

(2000) 

DEC. 

(2000) 

HEL.VEL. 

v(km/s) 

ERROR 

Av(km/s) 

Quality flag 

Emission lines 

151 

ESOl 

04:54:17.57 

-10:18:00.56 

75611 

88 

1 

NO 

152 

ESOl 

04:54:17.93 

-10:18:04.53 

-2 

-2 

1 

NO 

153 

ES02 

04:54:18.17 

-10:10:30.22 

74556 

67 

1 

NO 

154 

ES02 

04:54:18.19 

-10:13:38.63 

72587 

75 

1 

NO 

155 

ESOl 

04:54:18.39 

-10:13:22.62 

31052 

113 

3 

NO 

156 

ES02 

04:54:18.43 

-10:16:50.73 

88515 

91 

1 

on 

157 

ES02 

04:54:18.56 

-10:13:39.25 

97693 

43 

1 

OII,H/3,OIIIa,b 

158 

ESOl 

04:54:18.87 

-10:18:12.17 

74546 

75 

1 

NO 

159 

ES02 

04:54:18.87 

-10:11:42.93 

72600 

121 

1 

on 

160 

ES02 

04:54:18.92 

-10:15:16.89 

74785 

62 

1 

NO 

161 

ES02 

04:54:19.08 

-10:16:35.20 

63328 

96 

2 

NO 

162 

ESOl 

04:54:19.23 

-10:15:19.14 

92338 

88 

2 

NO 

163 

ESOl 

04:54:19.32 

-10:16:45.17 

109151 

74 

1 

on 

164 

ESOl 

04:54:19.75 

-10:13:49.60 

-2 

-2 

1 

NO 

165 

ES02 

04:54:19.78 

-10:11:27.08 

75907 

87 

1 

NO 

166 

ES02 

04:54:20.06 

-10:13:28.71 

73571 

91 

1 

NO 

167 

ESOl 

04:54:20.53 

-10:10:22.20 

-2 

-2 

1 

NO 

168 

ESOl 

04:54:20.65 

-10:17:26.06 

99128 

89 

2 

NO 

169 

ES02 

04:54:20.92 

-10:13:55.87 

-2 

-2 

1 

NO 

170 

ESOl 

04:54:20.96 

-10:16:44.77 

78326 

146 

1 

on,n/3 

171 

ES02 

04:54:21.12 

-10:09:39.91 

-2 

-2 

1 

NO 

172 

ESOl 

04:54:21.29 

-10:12:07.41 

101738 

64 

2 

NO 

173 

ESOl 

04:54:21.30 

-10:14:46.26 

55236 

86 

2 

NO 

174 

ES02 

04:54:21.63 

-10:12:02.65 

-1 

-1 

4 

NO 

175 

ESOl 

04:54:21.65 

-10:16:41.02 

73571 

64 

1 

NO 

176 

ESOl 

04:54:22.20 

-10:14:50.23 

-2 

-2 

1 

NO 

177 

ES02 

04:54:22.32 

-10:16:25.27 

73694 

51 

1 

NO 

178 

ESOl 

04:54:22.35 

-10:17:13.40 

72593 

107 

1 

H/3 

179 

ESOl 

04:54:22.64 

-10:16:40.67 

79233 

125 

1 

NO 

180 

ES02 

04:54:22.70 

-10:15:13.01 

54106 

71 

1 

n/?,onia,b 

181 

ES02 

04:54:23.18 

-10:17:58.06 

72568 

93 

1 

on,onia,b,n/? 

182 

ESOl 

04:54:23.19 

-10:12:38.32 

-2 

-2 

1 

NO 

183 

ES02 

04:54:23.25 

-10:11:33.53 

88370 

59 

1 

NO 

184 

ESOl 

04:54:23.73 

-10:13:05.92 

109876 

144 

2 

NO 

185 

ESOl 

04:54:23.73 

-10:17:10.67 

78122 

100 

1 

H/3 

186 

ES02 

04:54:23.93 

-10:11:26.49 

88524 

50 

1 

NO 

187 

ESOl 

04:54:24.16 

-10:14:00.74 

-2 

-2 

3 

NO 

188 

ESOl 

04:54:24.55 

-10:16:40.88 

73983 

48 

1 

NO 

189 

ESOl 

04:54:24.60 

-10:14:12.10 

106221 

88 

1 

on,n/3,onia,b 

190 

ESOl 

04:54:24.77 

-10:14:14.81 

75310 

84 

3 

NO 

191 

ESOl 

04:54:24.91 

-10:16:58.01 

72929 

69 

1 

NO 



